Levy flights in quantum transport in quasi-ballistic wires. 
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Conductance fluctuations, localization and statistics of Lyapunov exponents are studied numerically 
in pure metallic wires with rough boundaries (quasi-ballistic wires). We find that the correlation 
energy of conductance fluctuations scales anomalously with the sample dimensions, indicating the 
role of Levy flights. Application of a magnetic field defiects the Levy fiights which reduces the 
localization length. This defiection also breaks the geometrical fiux cancellation and restores the 
usual Aharonov-Bohm type magneto-conductance fiuctuations. 
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During recent years theoretical interest in mesoscopic 
systems has started to divert from studies of universal 
properties of ergodic systems, such as disordered metals, 
towards structures which are non-ergodic |l| . Ergodicity 
of diffusive conductors is provided by effective mixing of 
the phase space corresponding to the classical counter- 
part of quantum mechanics - due to quantum scattering 
at a short-range random potential. In non-ergodic sys- 
tems, homogeneous mixing of classical phase space does 
not occur before the quantum state is formed. Although 
issues of non-ergodicity are usually addressed in a context 
of spectral and wave function statistics in closed systems, 
such as partly chaotic billiards 0-0], one may extend 
the notion of non-crgodicity to open systems and study 
its implications in quantum transport effects such as lo- 
calization and conductance fluctuations. We distinguish 
below a class of disordered objects where a large part of 
classical phase space is not mixed p] since the particle 
motion through it is not diffusive. These systems are 
pure metallic wires with rough edges corrugated at the 
length scale of the electron wavelength (quasi-ballistic 
wires) j^-p^, where the existence of anomalously long 
grazing ballistic paths is known to lead to noticeable ef- 
fects even in classical transport properties |1J]. 

In the theory of dynamically generated random walks, 
ballistic paths with lengths very much exceeding the av- 
erage have the name Levy flights [ [l5[[lq| and when the 
latter occur with a relatively high probability they de- 
termine a superdiffusion character of random motion. In 
diffusive wires with static short-range bulk impurities, 
the probability of finding a ballistic path rj longer than 
the mean free path, / is exponentially small, ^ exp(~77/^). 
This leads to the normal diffusion relation between the 
variance of the lengths of random walk paths and the 
system size and gives the 'ergodic' form for the corre- 
lation energy of conductance fluctuations, E^ ^ hD/L^ 
m . A natural example of a classically superdiffusive sys- 
tem with static disorder is a pure metallic wire with cor- 
rugated edges, where Levy flights are just those graz- 
ing ballistic trajectories which cross it at a small angle, 



Q <C 1, to the wire axis. In wires with only edge dis- 
order, ballistic segments rj — W/ sina much longer than 
the sample width appear with the probability ^ 77"^ for 
T] 3> W. Below we report the results of numerical simu- 
lations of quantum transport in quasi-ballistic wires with 
edges which are rough on an atomic scale. We analyze 
the effects of superdiffusion on the correlation properties 
of mesoscopic conductance fluctuations, 6G{e, B) and the 
effect of a magnetic field on the localization length, L^- 
To anticipate a little, the correlation energy of conduc- 
tance fluctuations we found scales anomalously with the 
system size as E^ — ^hWvpL^^ ln( -^ ^^/ )- Although geo- 
metrical flux cancellation Iq-pQ] hinders the conventional 
crossover in the localization effects between the orthog- 
onal and unitary symmetry classes |0,18|, we observe 
the usual Aharonov-Bohm magnetoconductance fluctua- 
tions due to the deflection of Levy flights by a magnetic 
field (which transforms part of them into skipping or- 
bits). This deflection breaks the flux cancellation rules 
and also determines the magnetic field dependence of the 
localization length which is different to that obtained in 
ergodic (with bulk impurities) disordered systems. 

The results presented below are based upon numerical 
solution of a two-dimensional Anderson Hamiltonian on 
a square lattice, H ^ J2i N)<^i(«l - ^S{ii> K)0'l' where 
(ij) denotes nearest neighbour sites i and j. The struc- 
ture considered consists of two ideal leads of width W 
attached to a scattering region W sites wide and L sites 
long (all lengths are in units of the lattice constant a) . In 
the absence of a magnetic field, the off-diagonal matrix 
elements V — 1 determine the width of the energy band. 
Within the leads and in the centre of the scattering re- 
gion (when modeling clean wires with rough boundaries), 
the diagonal matrix elements were e^ = eo (with ep > 1, 
which keeps the Fermi level away from the van Hove sin- 
gularity in the band center). For simulating bulk disor- 
der, the ei in the scattering region are taken uniformly 
from the interval — C//2 < e^ — eg < U/2 where U is the 
disorder strength. For sites on the boundary, e^ = eo + es 



where cb — 10''. The rough structure of the boundary 
was generated by having an equal probability of either 0, 
1 or 2 sites at each edge with on-site potential eo + e^ 
[ [19[ . To simulate the effect of a magnetic field, we in- 
corporate a Peierls phase factor into V in the scattering 
region. Numerical decimation techniques pO| were used 
to compute the Greens function for each given realization 
of disorder. From these the transmission matrix tt^ was 
found and diagonalized. Its eigenvalues, r„, were ana- 
lyzed statistically and used to compute the conductance 
G = ^„ T„, which we measure in units of ^. Note that 
numerical results presented here describe a model system 
without any decoherence {i.e., without any inelasticity). 

The effect of Levy flights can be easily identified in 
our numerical results of quantum transport in pure wires 
with edge disorder. Rather than consider the statistics of 
the TnS it is more useful to introduce the parametriza- 
tion T„ = 1/cosh {L/S,n) and consider the distribution, 
P(^~^), of inverse localization lengths or Lyapunov ex- 
ponents. This is shown in Fig. 1 (a) for 4 series of quasi- 
ballistic wire samples with W — lb (giving a mean width 
of 13) and lengths L = 52 (A), 104 (B), 208 (C) and 
416 (D), and a series of samples (U) with bulk 'defects'. 
As pointed out by Tesanovich et al 0] and verified nu- 
merically in Ref. |ll[], the length of Levy fiights, r/max, 
is limited in quantum systems. This is because uncer- 
tainty in the transverse momentum, Sk± ~ ^/W, in a 
wire with a finite width sets a quantum limit to the an- 
gles W/r] ^ a > 6a ^ 5k±/kp ^^ X/W which can be 
assigned to a classically defined ballistic segment [nOJ. 
This sets the cut-off ?7max = W^/Xp, and a finite local- 
ization length Lc ~ ■^^x~~ ■ Samples from the series A and 
B meet the criterion L < Lc and the distribution of the 
Lyapunov exponents obtained in them has a pronounced 
peak at small ^~^ corresponding to eigenvalues T„ « 1. 
This should be compared to the plateaux-like [2^ dis- 
tribution P(^~^) reproduced using the same numerical 
procedure in a sample from the series U. Note a finite 
width of the balhstic peak in P{^~^) and that P(0) = 0. 
The enhanced density of small £,~^ can be identified even 
in samples from the series C and D with L k, L^ where 
P(^~^) starts to show a periodic modulation specific to 
the localized regime where the spectrum of Lyapunov ex- 
ponents tends to crystallize [ p^|j2^ . 

The distribution of the eigenvalues of the transmission 
matrix results in a finite-width distribution of conduc- 
tances and, hence, conductance fiuctuations. The statis- 
tics of conductance fluctuations, P(G'), for wires with 
edge disorder is shown in Fig. 1(b) (for the case of zero 
magnetic field) for the series A, VF = 15, L = 52 < Lc ~ 
ttW'^IXf ^ 200, (G) w 2.0 (the distribution is almost 
Gaussian) and G,W = \b, L = 208 ~ Lc, (G) « 0.71. 
The distribution function P{G) is the result of the anal- 
ysis of various realizations of disorder. We shall denote 
such averaging over realizations by (....). When calcu- 



lating correlation functions, we also add averaging over 
energy. The sample length dependence of the correlation 
energy, Ec{L), of the fluctuation pattern found from the 
half-height of the correlation function {SG{e)6G{e + 6e)) 
is shown in Fig. 1(c). For comparison, Ec{L) ^ L~^ is 
also shown for bulk-disordered samples (U) , which serves 
us as a reference case. 

In wires with edge disorder ballistic segments 77 = 
W/ sin a much longer than the sample width appear with 
the probability n-^\^^\-^ -- r]'^ for L > 77 > IF > 
Xf and this gives rise to a difference of the correlation en- 
ergy sample length dependence from the aforementioned 
ergodic law. The correlation energy is determined by 
the lowest eigenvalue in the relaxation spectrum of the 
diffusion operator. In the anomalous case, the diffusion 
process is described by the integral equation 



ft ,\ 1 /-w PU-^l + {x-y)^;y,x' 
P{--X,x']^- —^ :-j:: ^-dy 

with initial conditions P{t ~ 0;x,x') ~ d{x — x'). The 
correlation function P{t/T;x,x') describes the probabil- 
ity density to find a particle at the diffusively scattering 
edge of a sample symmetrized over both wire edges, and 
it is represented in the coordinates along the wire axis 
normalized by the wire width. The time is normalized 
by the ballistic time, r = W/vp. The solution of this 
equation in a long system can be found by taking the 
Fourier transform, 
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To study the low- lying spectrum of P(w,g), we expand 
n over u! ^ 1, and find that, 

n « (1 - qK^iq)) - — + ^ + - 

where Kn are Bessel functions. In the limit oi q -^ I this 
gives 
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For a sample of length L, we take q^nin = ttW/L and get 
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where 7 is Euler's constant. The anomalous dependence 
in this analytical result describes a faster escape of a par- 
ticle from a quasi-ballistic wire, as compared to a bulk- 
diffusive one, and it is represented by the solid curve in 
Fig. 1(c). Note that to relate the above semiclassical 
analysis to purely quantum numerical simulations, we 
use an effective sample width (W) — Xp. In addition. 



it is possible to calculate the variance of conductance 
fluctuations in a quasi-ballistic case, (SC^y,)- This has 
been done by evaluating the probabilities of two paths 
to encounter each other at two scatterers, which uses 
the result of Eq. (H). Such a calculation shows that 
('^^qb)/(^^erg) — |t ~ 0.66, whcrcas the numerical re- 
sult obtained for a system with W — 15 and L = 104 is 
{6Gl^)/{SGl^)^0.72. 

Another feature of quasi-ballistic wires is flux cancella- 
tion [p|-pO| , which is an exact geometrical fact in ballistic 
systems with parallel edges. Because of this cancellation, 
no Aharonov-Bohm flux can be encircled by loops com- 
posed of a sequence of ballistic flights between sample 
boundaries. In the metallic regime, L < L^, the curving 
of an electronic trajectory by a magnetic fleld reduces the 
length of the longest ballistic paths, so that the conduc- 
tance (G) declines when the transverse deflection of the 
Levy flights, Sx± fa ^{evpB/mc) (L/vp) ^ BL'^Xp/tpo 
becomes larger than the wire width, Sx± > W: 

at B > Bd.fi^ cI)qW/{XfL^) ^ {LjL){WLy^(bo. (3) 

This effect is present both in the averaged (over disor- 
der) conductance for samples from series A and B shown 
in Fig. 2(a) and in the suppression of a peak at ^^^ 
in P(f^^) in Fig. 2(b) for samples from the same se- 
ries, as in Fig. 1(a). The conductance (G) decreases up 
to the field where the electron cyclotron radius Re be- 
comes commensurable with the wire width W. At such 
a field, Scorn = 0.55-jpp- Q, the magnetoconductance 
has a minimum followed by a rise and formation of the 
quantum Hall effect plateaux after which any electron 
cyclotron orbit either decouples from the edges, or skips 
independently along each of them pj . 

The shortening of the length of Levy flights reduces the 
localization length in a wire. In samples from the series 
C and D (also shown in Fig. 2(b)), localization appears 
to be stronger in a field than at i? = 0: The crystal- 
lization in the spectrum of ^~^ prevails, which refiects 
the fact, supported by the analysis of longer wires (up to 
L — 1040) , that the localization length is shortened by a 
field. This tendency, illustrated in Fig. 2(c) ( soHd curve), 
is opposite to the behavior of Lc{B) ^^^^ in diffusive 
samples from the U-series (dashed curve). The relevant 
scale of magnetic fields at which the tendency of the lo- 
calization length to shorten in quasi-ballistic wires starts, 
(J)oXf/W^, can be obtained from Bdcfl by replacing the 
sample length L with Lc, which formally gives a similar 
scale as a field which would provide the Aharonov-Bohm 
type of a crossover in Lc between two symmetry classes, 
so that the latter is hindered by the deflection effect. 

The same deflection of Levy flights at B > Bdcfl also 
breaks the geometrical flux cancellation and restores the 
ability of a pair of electron paths to encircle a magnetic 
flux ~ LWB. At this fleld the longest ballistic flights, 
which still dominate transport, transform into skipping 



orbits so that the electron path inside a sample acquires 
segments which do not cross the wire: They start and end 
on the same edge. On the one hand, the positive magne- 
toconductance, which could result from the suppression 
of the weak localization correction developed at the field 
scale oi B ^ -Bdefl, is mostly hindered by a larger negative 
classical effect. On the other hand, the break-down of ge- 
ometrical flux cancellation manifests itself in pronounced 
magnetoconductance fluctuations 6G{B) = G{B) — (G) 
illustrated by an example in Fig. 3(a) (solid line), which 
shows the conductance G{B) of a single sample, along 
with results from a sample from the U-series (dashed 
line - this has been shifted by 1 vertically for clarity). 
The variance of fluctuations is shown in Fig. 3(b). From 
analyzing the auto-correlation function as a function of 
magnetic fleld, we deduce the sample length dependence 
of the correlation magnetic fleld Be in quasi-ballistic sam- 
ples. It corresponds to the magnetic fleld flux through 
the sample area equal to 2.5(/)o, in comparison with 1.500 
that we get in the bulk-disordered case. This seems to 
explain an earlier experimental observation ||24[| . 
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FIG. 1. (a) Distribution function P{L/^) for 
4 quasi-ballistic samples with W =15 and L =52(A), 104(B), 
208(C) and 416(D) and for one wire (U) with on-site bulk 
disorder U = 2.0 (nominal mean free path I ~ 8.5, H^ = 15 
and L = 52). (b) shows P{G) for the two quasi-ballistic struc- 
tures A and C. (c) shows Ec{L) for quasi-ballistic wires with 
£0 £ [1-5, 1.7] (circles) and eo G [1-0, 1.2] (diamonds) and for 
a wire from the U series (squares) . Also shown is the analyti- 
cal result (eq. (2) - solid curve) with hvf = 12.3 and W — 9.3 
[15]. The dashed curve is the ergodic law Ec ^ hD/L^ with 
hD = 62.5. 




2.0 4.0 6.0 8.0 10.0 




I ' ' ' ' ' ' ' 1 

0.000 0.005 0.010 0.015 0.020 

B (^o/a') 

FIG. 2. (a) The ensemble averaged conductance as a func- 
tion of a magnetic field for wires from the A and B series, 
(b) The same curves (A-D) as for Fig. 1(b) but at finite 
magnetic field B = 0.0200 /a^- (c) the localisation length for 
quasi-ballistic (solid curve) and disordered (dashed curve - 
bulk mean free path I ~ 8.5) wires. 
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3(a) Conductance for one quasi-ballistic wire (solid 
curve) from the B series and one disordered wire (dashed 
curve - mean free path / ss 16). 3(b) Variance var(G(B)) 
for quasi-ballistic wires from the B series (solid curve) 
and disordered wires (dashed curve - mean free path 
?«16) 



